Landau level mixing in the v — 5/2 fractional quantum Hall state 
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The v = 5/2 fractional quantum Hall state is studied numerically, directly including the effects of 
electron scattering between neighboring Landau levels. Significant reduction of the excitation gap 
caused by the LL mixing explains the discrepancy between earlier calculations and experiments. On 
the other hand, LL mixing also considerably reduces overlaps with the Moore-Read wavefunction, 
raising a question of the actual realization of nonabelian quasiparticles in present experiments. 
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I. INTRODUCTION 

Discovery 1 of the fractional quantum Hall (FQH) effect 
in a half-filled first excited Landau level (LLi) demon- 
strated the possibility for incompressible quantum liquid 
(IQL) states outside of the "standard" hierarchy 3 , de- 
scribed below. Laughlin states^ occur at the fractional 
LL fillings v = 2irg\ 2 = (2p + 1) _1 (where g is electron 
concentration, A is the magnetic length, and p is an inte- 
ger). The simple form of Laughlin many-body wavefunc- 
tions suggested that their low energy could be understood 
in terms of the avoidance of pair states with the smallest 
relative pair angular momentum (and largest repulsion); 4 
The quasiparticles (QPs) of the Laughlin states can (un- 
der some conditions^) form "daughter" IQL states with 
their own QPs, giving rise to an entire IQL hierarchy. 3 
The most stable IQL states, occurring at v = s(2ps±l) _1 
(s being another integer), appear naturally in Jain's com- 
posite fermion (CF) model 6 involving the concepts of flux 
attachment and an effective magnetic field. All Laughlin 
and Jain states are characterized by odd-denominator fill- 
ing fractions v and fractional QP charge q = e(2ps±l) _1 . 

For the even-denominator IQL state observed &tu — | , 
Moore and Read (MR) proposed 7 a different, paired 
wavefunction, and predicted that its QPs obeyed non- 
abelian statistics. The MR state has been studied in 
great detail^ and interpreted by two complementary 9 
pictures: as a Laughlin state of tightly bound electron 
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or a superfluid of weakly bound CF pairs 
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pair; 

The first numerical calculations for interacting elec- 
trons in a partially filled LLi were carried out by Morf. 14 
They seemed to confirm that a half-filled LLi has a 
spin-polarized incompressible ground state accurately de- 
scribed by the MR wavefunction. However, the subse- 
quent experimentsi£*i£*i2*i& revealed minute excitation 
gaps A ~ O.i — 0.45 K in the real v = I states, up to 
20 x smaller than predicted from numerics. The fact that 
only a small part of this discrepancy could be attributed 
to the finite width of the quasi-2D electron layer or to 
the weak disorder posed a challenge for the theoretical 
models. It also raised a fundamental question of whether 
the actual, experimentally realized u = | states are in- 
deed adequately described by the MR wavefunction. It 
has become very important in the context of topologi- 



cal quantum computation, whose recent proposals^^SiSi 
take advantage of the nonabelian QP statistics. This 
fundamental question is the main subject of our present 
paper. 

We directly include LL mixing in exact diagonalization 
by adding states containing a single cyclotron excitation 
to the Hilbert space of a partially filled spin-polarized LL. 
This can be regarded as a first-order approximation in 
the Coulomb-to-cyclotron energy ratio (3. We apply this 
procedure to the v = |, |, and | states and evaluate the 
excitation energy gaps A and (in the last case) the over- 
lap £ with the MR wavefunction. At (3 ~ 1 (relevant for 
the experiments at v = f), where our first-order results 
are only qualitative, we find a strong reduction in both 
A and £. Unfortunately, the calculations including exci- 
tations to higher order in (3 are beyond our capabilities. 
Our conclusion concerning the gap reduction gives quali- 
tative support to the work of Morf and d'Ambrumenil, 22 
who however included LL mixing in a different wayjS 3 . 
using a screened interaction that is strictly valid only 
for v ^> 1. Furthermore, our prediction of strong de- 
viation from the MR wavefunction supports the recent 
proposals— for new experiments aimed at determining 
directly the QP statistics at v = | (rather than merely 
at determining the incompressibility). 



II. MODEL 

We use Haldane's spherical geometry? 3 " convenient for 
the exact study of liquid states with short-range correla- 
tions. On a sphere of radius R, the normal magnetic 
field B is produced by a Dirac monopole of strength 
2Q = 4irR 2 B/<j> , defined here in the units of flux quan- 
tum </)q = hc/e. Using a magnetic length A = ^Jfic/eB, 
this can be rewritten as QX 2 = R 2 . The series of LLs 
labeled by n = 0, I, 2, ... are represented by shells of 
angular momentum I = Q + n and degeneracy g = 21 + 1. 
The cyclotron energy is nhuj c = heB / ' [ic (counted from 
the lowest LL), where /x is the effective mass. The or- 
bitals ipnm{S,4>) are called monopole harmonics. 

The A/"-electron Hamiltonian matrix is calculated in 
the configuration-interaction basis . . . , itf). Here, the 
composite indices i — [n, m, a] also include spin, and the 
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FIG. 1: (color online) (a) Energy spectrum (energy E as a 
function of angular momentum L) of N = 12 electrons in a 
LL shell of angular momentum I = 25/2, corresponding to 
an isolated, half-filled first excited LL. (b) Excitation gap A 
extracted from the spectra similar to (a) , plotted as a function 
of the electron number N; w is the width of the 2D electron 
layer and A is the magnetic length. 

expressions for two-body Coulomb matrix elements (also 
in layers of finite widths w) can be found for example in 
Ref. [25[ At a given LL filling (defined by Af and g, the 
basis states can be classified by AS — Tilu c J2k n k ~ £min, 
i.e., the total cyclotron energy measured from the lowest 
possible value £ m [ n allowed by the Pauli exclusion prin- 
ciple (e.g., £ m in = at v < 2, or £ m i n = (J\f — 2g)Twj c at 
2 < v < 4). Alternatively, the total number of cyclotron 
excitations K = A£/(hcu c ) can be defined. 

The exact numerical diagonalization in the Hilbert 
space restricted to K — means including Coulomb scat- 
tering within only one, partially filled LL, and neglecting 
the LL mixing. For v — | this reduces the AA-particle 
problem to N = J\f — 2g electrons confined to an isolated 
LLi. A typical numerical spectrum is shown in Fig.^a). 
The non-degenerate ground states with a gap generally 
appear in finite-size systems with even values of N (they 
are known to be paired^) at 21 = 2iV+l or (equivalent via 
the N —> g — N particle-hole conjugation) at 21 — 27V — 3, 
both extrapolating to N/g — > | for large N. As shown 
in Fig.^b), the excitation gap A rather weakly depends 
on N, allowing one to estimate the value A rs 0.02 e 2 /A 
for an infinite (planar) systemii&ii 

Unfortunately, this value is not confirmed by the ex- 
periments. The gaps measured from the thermal acti- 
vation of longitudinal conductance range from 0.001 to 
0.004 e 2 /A, depending on the electron mobility) 15 i 16 i 17 i 18 
with extrapolation to a disorder-free system not exceed- 
ing ~ 0.006 e 2 /A. As shown in Fig. ^b) for w — A (i.e., 
w = 11 A nm at B = 5 T), this discrepancy cannot be 
explained by a finite width of the electron layer. 

An obvious advantage of the K — approximation is 
that calculations can be done for sufficiently large values 
of N to eliminate finite-size errors. It could be trivially 
justified by a small ratio of Coulomb and cyclotron en- 
ergies, (3 = (e 2 X~ 1 )/(huj c ) — A/cib (with the Bohr radius 
a-B = h 2 /fie 2 ). However, (3 > 1 at the fields B ~ 5 T 
typically used in FQH experiments at v = I. 



FIG. 2: (color online) (a,b) Comparison of inter-LL excita- 
tions with up to one unit of cyclotron energy in the v = | 
state and in Laughlin/Jain liquids at v < 1. (c) Reduction 
of the excitation gap of different quantum Hall states due to 
LL mixing, plotted as a function of the cyclotron-to-Coulomb 
energy ratio /3 . A and Ao are gaps calculated for N = 8 
and 21 = 16, 21, and 17 (for v = §, |, and §). The full dots 
correspond to the electron concentration g — 2.3 x 10 11 cm -2 , 
(d) Gap reduction 1 — A/Ao as a function of parameter /3* 
in which Coulomb energy q 2 Xq 1 involves different fractional 
charge quanta q (indicated) appropriate for different liquids. 



To include LL mixing, we expand the Hilbert space by 
adding the K = 1 states, i.e., allow excitation of up to 
one electron to a higher LL. This can be regarded as a 
kind of first-order perturbation scheme in (3. However, 
note that such coupling of the entire K = and K = 1 
subspaces is more than a first-order perturbation on the 
specific K — eigenstates. The basis is shown schemati- 
cally in Fig.[5Ja) and (b). At v = |, the underlying, filled 
lowest LL gives rise to more types of excitations than at 
v < 1. Generally, the inter-LL excitations can be decom- 
posed into addition of an electron or a hole to a specific 
LL in the presence of a correlated state of the initial N 
electrons. There are two distinct cases, depending on 
the target LL: (i) Addition to (or removal from) the LL 
occupied by the incompressible liquid. This causes cre- 
ation of several fractionally charged QPs. At v = i or | 
they are well known Laughlin or Jain QPs, with a sim- 
ple and intuitive picture in the CF model. However, the 
nature and dynamics of the QPs at v = | are not nearly 
as well understood, (ii) Addition to (or removal from) a 
different LL. This makes the added/removed electron dis- 
tinguishable from the correlated electrons. This problem 
resembles coupling of a Laughlin liquid to a (positive) va- 
lence hole 2 ^ or a (negative) trioni^l However, coupling of 
the v = | state to a foreign charge is far less understood. 

Since all three types of K — 1 excitations must be in- 
cluded in the calculation on the same footing, even this 
limited account of LL mixing boosts the space dimension 
from - 10 3 to ~ 4 x 10 5 for N = 8 at 21 = 17. This pre- 
cludes similar calculations for larger systems or further 
inclusion of the K > 1 excitations. On one hand, this 
makes the present results somewhat susceptible to finite- 
size errors (although Fig.^b) may suggest that N = 8, 
i.e., four pairs, is already a representative system) and 
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accurate only up to the first-order perturbation in /?. On 
the other hand, a much larger number of K = 1 excita- 
tions at v > 2 than at v < 1 suggests that the effects of 
LL mixing should be more important at v = | than in 
Laughlin or Jain liquids of the lowest LL. 

III. RESULTS AND DISCUSSION 

We focus on two features of the v = § state: the ex- 
citation gap and the overlap with the MR wavefunction. 
The results of calculations of the gap A for N = 8 and 
K < 1 in a 2D layer of zero width are shown in Fig. Et c )- 
As argued above and anticipated from experiments, the 
gap reduction A/ An (where Ao is the result for K = 0) is 
noticeably greater at v = | than at both v = | or |. On 
the other hand, the fact that the gap at v = ^ is reduced 
more than at v = | can be related to the smaller QP 
charge in the latter case (q = e/5 versus e/3). Since the 
low-energy response of a liquid involves formation and 
interaction of the QPs, the harsh /3< 1 criterion for the 
accuracy of the isolated-LL approximation may be re- 
laxed to (3* < 1. Here, (3* = (q 2 \~ x )/(?iuj c ) = {q/ef/ 2 (3 
involves the Coulomb energy scale of the QPs. Indeed, 
when the gaps in the weak-perturbation regime are plot- 
ted as a function of (5* as in Fig.^d), the data for v = i 
and | falls close to the same line, 1 — A/A ~ 3/3*. Tak- 
ing q — e/4 for the v = | state results in a much (about 
3x) steeper curve. This indicates that the response of 
the v = | state to the perturbation associated with the 
LL mixing at a finite (3 is (due to a richer inter-LL exci- 
tation spectrum) relatively stronger than the response of 
Laughlin or Jain states in the lowest LL. 

Furthermore, if the experiments on all three electron 
liquids were to be carried out at similar concentrations 
(corresponding to a maximum mobility), the difference 
between them will be additionally magnified by a differ- 
ence in f3 corresponding to different v. For example, for 
g = 2.3 x 10 11 cm~ 2 we obtained gap reduction of 8.5%, 
2.5%, and 35% at v = A, |, and |, respectively. 

In the above discussion we have established the follow- 
ing, (i) Realistic estimates of the excitation gap at v = | 
must include the LL mixing, whose effect at this filling 
is much stronger than for Laughlin or Jain states in the 
lowest LL. (ii) The gap reduction caused by LL mixing is 
already significant in the K < 1 approximation, (iii) It 
is plausible that the full account of the LL mixing might 
reconcile experimental results in the limit of vanishing 
disorder with the numerics. Unfortunately, calculations 
for K > 1 and N > 8 are beyond our present capabilities. 

The dependence of gap A on the magnetic field B and 
finite layer width w are displayed in Fig. |3J Finite width 
was introduced to the model by the calculation of two- 
body Coulomb matrix elements using 3D wavefunctions 
x{z)^ nm (9,(j)) with x( z ) cos(7rz/u;) for the normal 
direction. Remarkably, A(w) is nonmonotonic, with a 
maximum between w = 5 and 10 nm, depending on B. 
While the reduced gaps are still about twice larger than 



thin lines: calculation in isolated n=1 LL 
thick lines: including inter-LL excitations 
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FIG. 3: (color online) Dependence of the excitation gaps A 
calculated at v = | (for N = 8) with and without LL mixing, 
on the 2D layer width w (a) and the magnetic field B (b). 

the experimental values, the inclusion of even only K = 1 
excitations clearly improves the model. This suggests LL 
mixing as the main reason for the earlier A-discrepancy. 

Let us now turn to the question of equivalence of the 
v — | state realized in experiment and the model MR 
wavefunction. There is a subtle difference between the 
half-filled state and odd-denominator liquids like v = i 
or |. In the latter states, it is not merely incompress- 
ibility but also the form of correlations and many-body 
wavefunction that are robust against the variation of ma- 
terial, w, or B — as long as the interaction pseudopoten- 
tial is sufficiently strong (superharmonic) at short range, 
weak compared to the cyclotron energy, and strong com- 
pared to disorder. In contrast, the half-filled state re- 
mains incompressible for a wide class of electron-electron 
pseudopotentials, but the exact form of the wavefunction 
strongly depends on their details. This limits information 
that can be inferred about the nature of the state from 
the observation of its incompressibility. 

The underlying reason is the competition of at least 
two distinct wavefunctions sharing the same symmetry: 
the MR state, which can be defined as an exact ground 
state of a short-range three-body repulsion^ and a clus- 
tered state characterized^ by the maximum avoidance 
of the next to the lowest value of the relative pair an- 
gular momentum, TZ = 3. The MR state is anticipated 
for the pair repulsion which is nearly harmonic at short 
range (i.e., with the pseudopotential decreasing linearly 
through 1Z = 1, 3, and 5), such as in the LLi. How- 
ever, the overlaps of the actual Coulomb eigenstates ob- 
tained from finite-size numerics (in the K = approxi- 
mation) with the MR state are sensitive to the interaction 
parameters^ and surface curvature^i raising the ques- 
tion of whether the v — | FQH state and the MR model 
state are indeed (qualitatively) equivalent. 

How does LL mixing affect this problem? In the top 
frames of Fig.^we plot the squared projection \Vk=o^\ 2 
onto the LLi (i.e., onto the K = subspace) of the same 
ground states ^ whose gaps are shown in Fig. It de- 
pends on w and B, but it is always significantly higher 
than A/ An or the (not shown) squared overlap with the 
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thin lines: calculation in isolated n=1 LL 
thick lines: including inter-LL excitations 
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FIG. 4: (color online) Dependence of the squared projection 
onto the LLi (top) and the squared overlap with the Moore- 
Read wavefunction (bottom) of the u = | ground states cal- 
culated for N = 8 with and without LL mixing, on the 2D 
layer width w (left) and the magnetic field B (right). 



K = ground state, in consequence of the coupling 
between intra- and inter-LL excitations in the K < 1 
space. In the bottom frames we plot the squared overlap 
£ 2 = |(MR|*)| 2 with the MR state (more precisely, with 
its particle-hole conjugate at 21 = 2N+1). The small val- 
ues for the K — calculation (here, 0.75 to 0.80) may be 



to some extent an artifact of spherical geometry 11 How- 
ever, a significant drop caused by the LL mixing (e.g., 
from 0.78 to 0.43 for w = 20 nm and B = 5 T) suggests 
that the MR wavefunction may not be a very realistic 
description of the v = § state in this range of parame- 
ters. This ambiguity and the difficulty with more realistic 
calculations make further experiments^ irreplaceable. 



IV. CONCLUSION 



We studied the effect of LL mixing on the v = | FQH 



state. This effect is more pronounced than in the incom- 
pressible liquids of the lowest LL such as v — | or | 
(even at the same magnetic field) due to the richer inter- 
LL excitation spectrum. The resulting reduction of the 
excitation gap explains the troubling disparity between 
previous numerics (without LL mixing) and the experi- 
ments. This prediction agrees qualitatively with Ref. 113 
which included LL mixing in a different way. Finally, 
the LL mixing significantly lowers overlaps with the MR 
wavefunction. This amplifies the need for the recently 
proposed^* experiments, designed to probe directly the 
nonabelian statistics of the QPs at v = §. 
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